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PREFACE 


This  report  has  been  prepared  to  assist  AWJ  personne?  in  estimating  the 
climatological  probability  of  aircraft  icing  according  to  season,  altitude, 
and  area  within  the  Northern  Hemisphere.  Tne  procedure  and  data  used  to  de¬ 
velop  isopleth  charts  of  icing  probability  acre  described. 

These  charts  can  be  used  in  support  of  long-range  planning.  Information 
on  short-range  planning  and  day-to-day  forecasting  of  icing  co-  <tions  is 
contained  in  AWSM  2  Of,- 39,  '‘Forecasters '  Guide  on  Aircraft  Ic.  * 

The  assistance  of  Mr.  Clarence  B.  Elam,  Jr.  is  gratefull;  acknowledged. 


LAWRENCE  G.  KATZ 
Environmental  Analyst 
ETAC,  USA? 

Washington,  D.  C.  20333 
30  December  i960 
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CLIMATOLOGICAL  PROBABILITY  OP  AIRCRAFT  ICING 


SECTION  A  —  INTRODUCTION 


In  1957  the  National  Advisory  Committee  for  Aeronautics  [4]  presented  air¬ 
craft  icing  data  which  were  obtained  from  various  Air  Heather  Service  recon¬ 
naissance  flights.  The  data  were  collected  at  700  mb  and  500  mb.  These  data 
readily  lent  themselves  to  statistical  analysis  because  of  the  manner  in  which 
they  were  collected.  The  reconnaissance  aircraft  flew  every  day,  at  the  same 
levels  and  over  the  same  routes.  Furthermore,  conditions  conducive  to  icing 
were  neither  sought  out  nor  avoided.  An  analysis  of  these  data  were  provided 
by  Appletnan  [1]  in  1959.  The  results  of  this  analysis  are  presented  as  graphs 
which  depict  the  conditional  probability  of  aircraft  icing,  given  cloud 
amounts  >  6/10,  as  a  function  of  altitude  and  temperature. 

In  1964,  Ingram  and  Gulllor.1  presented  an  atlas  which  provided  estimates 
of  the  percent  frequency  of  potential  icing  conditions  at  various  levels.  Po¬ 
tential  icing  was  considered  to  be  the  product  of  two  factors  —  the  probabil¬ 
ity  of  temperatures  <  C'C  and  the  probability  of  having  a  cloud  amount  >  6/10. 
It  was  subsequently  realized  that  potential  Icing  grossly  overestimated  actual 
icing  frequency.  Hq  Air  Heather  Service  then  requested  that  the  Environmental 
Technical  Applications  Center  (ETAC)  provide  charts  of  actual  icing  probabili¬ 
ties  based  upon  the  best  data  available.  This  report  was  prepared  to  fulfill 
that  request. 


SECTION  B  —  DEFINITION  AND  ASSUMPTIONS 


Definition 

The  probability  of  aircraft  icing  at  a  specified  altitude  in  the  atmos¬ 
phere  is  the  chance  that  an  aircraft  at  that  specified  altitude  will  in  fact 
experience  icing.  Estimates  of  icing  probabilities  are  based  on  climatologi¬ 
cal  estimates  of  the  parameters  involved. 


1  Ingram,  D.  M.  and  Gullion,  J,  L.:  "Estimated  Frequencies  of  Potential  Icing 
Conditions  at  Specified  Altitudes, "  Technical  Report  182.  Kq  Air  Heather 
Service  (ETAC),  17  pp.,  September  1964  (superseded  by  AHSM  105-39.  22  August 
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Assumptions : 

a.  Probability  of  icing  with  cloud  amounts  less  than  6/10  is  negligible 
compared  to  that  with  cloud  amounts  ^  6/10.  The  condition  that  the  cloud 
amount  be  equal  to  or  greater  than  6/10  is  assumed  necessary  but  not  suffi¬ 
cient. 

b.  The  range  of  temperatures  necessary  for  icing  extends  from  -2*C  to 
~4o-C.  This  report  is  prepared  for  conventional  (reciprocating  engine)  air¬ 
craft  and  so  it  is  assumed  that  the  magnitude  of  dynamic  heating  (heating  due 
to  friction  and  adiabatic  compression  as  the  aircraft  penetrates  the  air)  is 
two  Centigrade  degrees.  Furthermore,  it  is  assumed  that  no  supercooled  water 
exists  at  temperatures  lower  than  -4o*C .  The  temperature  range  -2*C  tc  -4o*C 
is  assumed  necessary  but  not  sufficient. 

c.  Weighted  averages  of  means  and  variances  of  temperatures  at  constant- 
pressure  levels  can  be  used  to  estimate  the  same  statistics  at  constant-height 
levels. 

d.  Frequencies  of  temperatures  below  -2*C  at  any  specified  altitude  are 
assumed  to  be  statistically  Independent  of  frequencies  of  cloud  amounts 

6/10  at  the  same  altitude. 

e.  The  mean  cloud  amount  at  an  altitude  is  assumed  to  be  equal  to  the 
probability  of  cloud  amounts  ^  6/10  at  that  same  altitude.  This  assumption 
appears  Justified  when  one  studies  the  shape  of  the  cloud-amount  distribution 
at  any  specified  altitude.  In  most  cases  it  is  a  U-shaped  distribution; 
therefore,  only  a  small  error  is  introduced  with  this  assumption. 


SECTION  C  —  PROCEDURE 


Atmospheric  temperature  is  a  continuous  function.  The  probability  of  air¬ 
craft  icing  at  a  specified  altitude,  location,  and  month  is  defined  mathe¬ 
matically  by  the  following: 

-2*C 

(I)  P(I)  -  f  F[I(T)|0  ^  6/10)  P[C  ^  6/10]  F[T]d! 

-40*C 

where: 

Finnic  £  6/10]  -  conditional  probability  density  function  of  icing 

as  a  function  of  temperature  given  that  the  cloud 
amount  is  equal  to  or  greater  than  6/10. 
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Ptl(T)|C  i  6/10)  *  {J  -  T  S 


T  <  -2*C 


Equation  (2)  Implies  that  although  the  range  of  temperatures,  -2*C  to  -40*C, 
and  cloud  amounts  ^  6/10  are  necessary,  r,hey  in  themselves  are  not  sufficient. 
Therefore,  icing  will  not  o..cur  every  time  an  aircraft  is  in  clouds  of  amount 
>  6/10  and  the  temperature  is  between  -2*C  and  -40*C.  Other  factors  such  as 
aircraft  type,  cloud  type,  liquid-water  content,  size  distribution  of  cloud 
droplets,  etc.,  play  a  part  in  determining  those  cases  conducive  to  aircraft 
icing. 

?[C  >  6/10]  is  the  probability  that  the  cloud  amount  is  equal  to  or 
greater  than  6/10.  F[T]  is  the  probability  density  function  of  temperature . 

By  definition: 


(3) 


F[I(T)|C  ±  6/10] 


F[I(T),  C  >  6/10] 

ric  >  6/io] 


thus, 

-2*C 

(4)  P(I)  -  J  F[I(T),  C  ^  6/10]  F[T]dT 

-40*C 

For  computational  purposes,  Zquation  (4)  is  replaced  by  the  approximation: 

(5)  P(I)  -  £  Fll(T),  C  2.  6/10]^  F[T]1  AT^ 

1«1 

where: 

F[l(T),  C  2.  6/'.0]1  =  probability  of  icing  in  temperature  class  inter¬ 
val  &T^  with  a  cloud  amount  equal  to  or  greater 
than  6/10 

F[T]^  -  probability  of  temperatures  in  temperature  class 
interval  AT^ 

Equation  (5)  is  Justified  since: 
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(6) 

n  -2*C 

£  F[I(T),  C  Z  6/10^  P[T]X  AT1  -  f  F[l(T),  C  ^  6/10]  P[T]dT 
1-1  -40*C 


P(I) 


or. 


(7) 

n 

lim  Y  F(l(T),  C  6/10 ].  P[T]t  AT, 

AT.-*  0  /-* 

1  1-1 


-2*C 

-  J  F[ I(T),  C  >  6/10]  F[T]dT 
-40*C 


P(I) 


For  the  purposes  of  this  report,  a  temperature  Interval  (AT< )  of  two  Centigrade 
degrees  was  used.  Since  temperatures  between  -2*C  and  -4o*C  are  assumed  to  be 
statistically  Independent  of  the  frequency  of  cloud  amounts  ^  6/10  at  the  same 
altitude,  location,  and  month,  we  may  write  Equation  (5)  In  the  following 
manner: 

n 

(8)  P(I)  -  P(C  >  6/10]  y  F[I(T)|C  >  6/10^  F[T]1 

1-1 

From  AWSM  105-39  ([2]  page  33),  Figure  12  presents  a  graphical  representation 
of  F[ I(T) | C  >  6/10]  for  the  700-  and  500-mb  pressure  surfaces.  Figures  14 
and  15  (AWSM  105-39,  PP.  41-42)  are  an  outgrowth  of  Figure  12.  Presented  in 
these  two  figures  are: 
n 

y  FlI(T)|C  >  6/10 ]±  F[T]^ 

1-1 

for  the  700-  and  500-mb  pres 3 tire  surfaces.  An  Inherent  assumption  in  these 
two  figures  Is  that  the  temperatures  at  these  levels  are  normally  distributed, 
and  a3  such,  are  completely  defined  by  a  mean  and  standard  deviation. 

Graphs  similar  to  those  presented  In  Figures  12,  14,  and  15  (AWSM  105-39), 
were  produced  by  ETAC  for  the  850-mb  pressure  surface.  These  graphs  were  pre¬ 
pared  from  pll.-w  reports  of  icing,  cloud  cover  as  shown  on  weather-dcpictlon 
charts,  and  85G-mb  temperature  analyses .  The  pilot  reports  and  weather- 
depiction  charts  were  used  to  determine  the  percent  frequency  of  icing  and 
no-iclng  occurrences  In  ^  6/10  clouds.  The  850-mb  charts  provided  a  tempera¬ 
ture  for  each  of  the  cases.  Based  on  this  information  a  probability  density 
function  of  Icing  in  clouds  >  6/10  as  a  function  of  temperature  at  the  850-mb 
pressure  surface  was  derived.  Thus,  the  probability  of  icing,  P(I),  for  a 
specified  altitude,  location,  and  month,  can  be  determined  from  a  knowledge  of 
the  mean  temperature  and  corresponding  standard  deviation  and  probability  of 
cloud  amounts  >  6/10. 
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SECTION  D  —  DISCUSSION  OP  DATA 

Temperature  Data 

Temperature  data  from  approximately  300  Northern  Hemisphere  stations  were 
used  in  this  report.  Monthly  mean  temperatures  ard  corresponding  standard 
deviations  for  the  850-,  700- ,  and  500-mb  pressure  surfaces  were  available. 

At  each  altitude,  seasonal  values  for  the  mean  and  standard  deviation  were 
derived  as  follows: 

*  V3  l  * 

1  n-1  1 

where: 

T  «  mean  seasonal  temperature  at  altitude  i 
*i 

and 

T  «  mean  temperature  for  month  n  at  altitude  i 
ni 

The  winter  season  is  defined  as  December,  January,  and  February;  spring  as 
March,  April,  and  May;  summer  as  June,  July,  and  August;  and  autumn  as  Septem¬ 
ber,  October,  and  November. 

(10)  a 

where: 

a  =*  seasonal  standard  deviation  of  temperature  at  altitude  i 
31 

and 

an  °  standard  deviation  of  temperature  for  month  n  at  altitude  1 

The  months  used  for  each  season  are  the  same  as  those  specified  above.  The 
seasonal  mean  temperature  and  standard  deviation  at  5000  feet  was  assumed  to 
be  equal  to  those  statistics  at  the  850-mb  pressure  surface.  Similarly, 
10,000-foot  data  were  taken  directly  from  the  700-mb  data.  The  seasonal  mean 
temperature  for  the  15* 000-foot  level  was  computed  using  the  following 
weighted  average: 

(11)  Ta  “  °*4  Ts  +  °*6 

3 15, 000  s700  3 500 

where: 
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T  -  mean  seasonal  temperature  at  15,000  feet 

“15,000 

T  -  mean  seasonal  temperature  at  700  mb 

“700 

and 

T  «  mean  seasonal  temperature  at  500  mb 

®500 

The  seasonal  standard  deviation  of  temperature  for  15,000  feet  was  computed 
using  the  following  equation: 


(12) 


“15,000 


of  +  0.6 

“700 


where: 

a  x  seasonal  standard  deviation  of  temperature  at  15,000  feet 

®15,000 

a  >  seasonal  standard  deviation  of  temperature  at  700  mb 

*700 

and 

a  x  seasonal  standard  deviation  of  temperature  at  500  mb 

s500 

The  seasonal  mean  temperature  at  20,000  feet  was  estimated  by  the  following 
equation: 


(13)  T.  -  T_  +  AT 

8 20, 000  500 

where: 

T  x  mean  seasonal  temperature  at  20,000  feet 

S 20, 000 

T  x  me  am  seasonal  temperature  at  500  mb 

“500 

and 

AT  -  temperature  difference  between  500  mb  and  20,000  feet 

taken  from  the  ICAO  Standard  Atmosphere  *  -3.4*C 

The  seasonal  standard  deviation  of  temperature  at  20,000  feet  was  assumed 
equal  to  that  at  500  mb. 


Cloud  Data 

Cloud  data  for  the  5000-,  10,000-,  15. 000- ,  and  20,000-foot  altitudes 
were  computed  by  the  method  developed  in  AWSTR  186  [3] .  Mean  seasonal  cloud 
amount  data  for  106  stations  were  utilized  in  this  study.  The  method  appear¬ 
ing  in  AWSTR  186  gives  mean  cloud  amounts  for  specified  levels.  The  mean 
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cloud  amount  at  an  altitude  is  assumed  equal  to  the  probability  of  ^  6/10 
cloud  amount  at  that  same  altitude. 

SECTION  E  —  CHARTS 


Charts  of  tne  conditional  probability  of  icing,  given  a  cloud  amount 
>  6/10,  were  constructed  by  combining  the  mean  temperatures  and  corresponding 
standard  deviations  with  the  appropriate  graph  of  P[l(T)|C  >  6/10].  Grid 
values  were  then  taken  off  the  charts.  Charts  of  P[C  >  6/10]  were  then  drawn 
and  grid  values  taken  off.  The  resulting  grid  maps  were  then  graphically 
multiplied.  The  resulting  chart  is  the  probability  of  aircraft  icing  at  a 
given  level. 

An  attempt  has  been  made  tc  ind' "ate  topographical  areas  of  greater  alti¬ 
tude  than  the  altitude  of  the  chart  on  which  they  appear.  These  appear  as 
shaded  sections  on  each  map. 

The  accompanying  charts  provide  estimates  of  the  probability  of  aircraft 
icing  at  the  5000-,  10.-000-,  15,000-,  and  20,000-foot  altitudes.  Charts  are 
presented  for  the  winter,  spring,  summer,  and  autumn  seasons. 
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FIGURE  2. 
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PROBABILITY  OP  ENCOUNTERING  ICING  CONDITIONS 

WINTER,  15,000  FEET 

FIGURE  3. 
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PROBABILITY  OF  ENCOUNTERING  ICING  CONDITIONS 

WINTER,  20,000  FEET 

FIGURE  4. 
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PROBABILITY  OF  ENCOUNTERING  ICING  CONDITIONS 


SPRIN6,  5,000  FEET 


FIGURE  3 
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FIGURE  7. 
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PROBABILITY  Of  ENCOUfc  TCRIN6  ICING  CONDITIONS 

SUMMER,  15,000  FEET 


FIGURE  II. 
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PROBABILITY  OF  ENCOUN1  cRING  ICING  CONDITIONS 

SUMMER,  20,000  FEET 

FIGURE  12. 
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PROBABILITY  Of  ENCOUNTERING  ICING  CONDITIONS 

AUTUMN,  10,000  FEET 

FIGURE  14 
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PROBABILITY  OF  ENCOUNTERING  ICING  CONDITIONS 


AJTUMN,  15,000  FEET 

FIGURE  IS. 


Technical  Report  194 


January  1967 


\ 


PROBABILITY  OF  ENCOUNTERING  ICING  CONDITIONS 

AUTUMN,  20,000  FEET 

FIGURE  1C. 
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